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The gluon and ghost propagators are computed using both quenched and dynamical configura- 
tions for the same lattice spacings. The Wilson fermions are simulated at several quark masses. 
Furthermore, the effect of the Gribov copies is evaluated for all sets of configurations. 
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1. Introduction and motivation 

It is well-known that the Faddeev-Popov (FP) procedure ^ can not be applied to the non- 
perturbative regime of Quantum Chromodynamics (QCD), due to the so-called Gribov copies pro- 
blem. The FP trick assumes a unique solution of the gauge fixing condition on each set of gauge 
related configurations. Gribov ^ has shown, for the Landau and Coulomb gauges, that there are 
more than one configuration satisfying the gauge condition on each gauge orbit. Furthermore, it is 
not possible to define a local continuous gauge fixing condition free of Gribov copies [|3[ |J. 

Many (quenched) lattice simulations have confirmed the existence of Gribov copies [j|], and 
the influence of this issue on Landau gauge gluon and ghost propagators have been studied (see, 
for example, [Q, |7|]). The simulations show that for the gluon propagator the influence of Gribov 
copies is null or small, whereas for the ghost propagator the effect is larger. 

In this paper we recover the computation of gluon and ghost propagators on dynamical con- 
figurations, including an investigation of the Gribov copies. Similar studies have been performed 
by other groups while ago ^ but, to the best of our knowledge, no investigation on the influence 
of the Gribov copies in dynamical simulations is found in the literature. 



2. Lattice setup 

In this paper we will consider 16 4 configurations. Three different sets of dynamical configu- 
rations for Wilson gauge action at jS = 5.6, using Wilson fermions, with different fc's, see table Q, 
were considered. For the dynamical ensembles, the lattice spacing have been computed using the 
Sommer scale ro = 0.5 fm. Data for aro, m q and m n have been taken from []|]. The j8 values for 



the corresponding quenched ensembles have been computed using [10] so that they have the same 



lattice spacing as the dynamical counterparts. Quenched simulations have been performed using 



MILC code [11], whereas dynamical simulations have been performed using Chroma software li- 
brary [12]; the MD integrator used in the HMC simulations has been tuned using the techniques 
described in [|j~3|]. 

Details on the lattice Landau gauge fixing and computation of the propagators can be found in 
[J7|]. The lattice data for the gluon propagator has been renormalized by requiring 

D R (q 2 )\ q 2 = ^ = ^=Z R D Lat (^ 2 ) (2.1) 

at = 3.5 GeV. In what concerns the ghost propagator, due to the very small statistical errors, 
we were unable to fit the high momenta in order to perform the renormalization. Given that we 





Quenched 






Dynamical Wilson 




p 


a(fm) 


# conf . 


K 


K mass (MeV) q mass (MeV) 


# conf. 


6.07 


0.0828 


60 


0.1575 


665 66 


71 


6.0845 


0.0809 


66 


0.1580 


485 34 


71 


6.106 


0.0782 


62 


0.15825 


385 22 


71 



Table 1: Lattice setup. 
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Renormalized Gluon Propagator - u. = 3.5 GeV 

a = 0.0828 fm; 1/a = 2.383 GeV 
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(a) Gluon propagator at a = 0.0828fm. 



Renormalized Gluon Dressing Function - u. = 3.5 GeV 

a = 0.0828 fm; 1/a = 2.383 GeV 
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(b) Gluon dressing function at a = 0.0828fm. 
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Renormalized Gluon Propagator - [l = 3.5 GeV 

a = 0.0809 fm; 1/a = 2.439 GeV 



Quench. 6 = 6.0845 1 
Dyn, K= 0.1580 | 
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(c) Gluon propagator at a = 0.0809fm. 



Renormalized Gluon Dressing Function - u. = 3.5 GeV 



a = 0.0809 fm; 1/a = 2.439 GeV 
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(d) Gluon dressing function at a = 0.0809fm. 
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Renormalized Gluon Propagator - u, = 3.5 GeV 

a = 0.0782 fm; 1/a = 2.523 GeV 



Quench. B = 6. 106 
Dyn, K= 0.15825 



q [GeV] 

(e) Gluon propagator at a = 0.0782fm. 
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(f) Gluon dressing function at a = 0.0782fm. 



Figure 1: Dynamical gluon propagators and dressing functions versus quenched counterparts. 
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m q |MeV| MJNL-V] 

(a) -D(O) vs. q mass. (b) D(0) vs. n mass. 



Figure 2: Zero momentum gluon propagator as a function of quark and pion masses. 

observed no difference in the UV between quenched and unquenched propagators, we just show 
the bare data. 

In order to evaluate the effect of Gribov copies, we have generated 100 copies per gauge 
configuration. We will consider, besides a random copy, the copies which have the largest and 
smallest value of the gauge fixing functional. 

3. Results 

In figure [l] the gluon propagator computed with dynamical fermions is compared with the 
quenched propagator. Figure [l] shows clearly that the effect of dynamical quarks is mainly seen in 
the infrared region. With the exception of the zero momentum, the propagator for the dynamical 
configurations is suppressed in the infrared. On the other hand, D(0) from the dynamical ensembles 
is suppressed for heavier quarks, and enhanced for lighter quarks. Furthermore, D t / V „(0) follows a 
linear behaviour as a function of the quark and pion masses - see figure |2| 

In what concerns the ghost propagator, see figure |3j, the effect of dynamical configurations is 
again more important at low momenta, with the dynamical quarks enhancing the propagator in the 
infrared region, if compared with quenched results. The data also shows that the enhancement goes 
with the quark mass, i.e. the larger the quark mass, the larger is the enhancement. 

The influence of the Gribov copies in the propagators is resumed in figures || and ^. Indeed, 
the pattern observed for the quenched and unquenched theories is similar: no or small effects for 
the gluon propagator; larger effects for the ghost propagator. Note that the Gribov effects in the 
ghost propagator are seen over a large range of momenta, with the propagator being suppressed for 
larger F maxima. 

4. Conclusions 

We have presented a first study of the effect of Gribov copies in Landau gauge gluon and 
ghost propagators computed from dynamical configurations. The results show that the effect is 



4 



Unquenching the Landau Gauge Lattice Propagators and the Gribov Problem 



Paulo J. Silva 



very similar for quenched and dynamical propagators. Furthermore, we have revisited the effects 
of dynamical quarks in gluon and ghost propagators. 
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Bare Ghost Propagator 

a = 0.0828 fin; 1/a = 2.383 GeV 



Quench, fi = 6.07 
Dyn, K = 0.1 575 



(a) Ghost propagator at a = 0.0828fm. 



Bare Ghost Propagator 

a = 0.0809 fm; 1/a = 2.439 GeV 



Quench, p = 6.0845 
Dyn, K = 0.1580 



(b) Ghost propagator at a = 0.0809fm. 



Bare Ghost Propagator 

a = 0.0782 fm; 1/a = 2.523 GeV 



Quench, (3 = 6.106 
Dyn,K = 0.15825 
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(c) Ghost propagator at a = 0.0782fm. 



Figure 3: Dynamical bare ghost propagators versus quenched counterparts. 
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Gluon Ratios 
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(a) Quenched gluon ratios at /3 = 6.106. 



Gluon Ratios 

K=0.15825 
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(b) Dynamical gluon ratios at K = 0. 15825. 



Figure 4: Gluon ratios. 




Ghost Ratios 

K=0.1575 
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(a) Quenched ghost ratios at = 6.07. 



(b) Dynamical ghost ratios at K = 0.1575. 



Figure 5: Ghost ratios. 
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